FINAL STATE INTERACTIONS IN B -> Dp and B -> D*vr DECAYS 
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We analyze final state strong interaction effects in B — ► and .B — ► -D*7r 
decays using the Regge model. We find that, due to the smallness of the 
contributions from the non-leading Regge trajectories (p , f, tv etc.), final state 
interaction phases are small if the Pomeron coupling to the charm quark is 
suppressed in comparison to lighter quarks. Our conclusion is that for B decays 
into states containing charm, final state interaction effects should play a minor 
role. 
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The problem of final state strong interactions in non-leptonic heavy meson decays 
has recently received considerable theoretical interest p], f| || . The relevance of 
the problem is related to the present and future programs for studying CP- violation 
in heavy hadrons, in particular B meson, decays. As well known ||, CP- violation 
in such systems might be observed by measuring an interference effect between two 
different amplitudes, and the relevant physical observable, i.e the CP-odd asymmetry, 
turns out to depend crucially on the strong interaction phase difference between the 
two amplitudes. 

Whereas the usual approach for charmed meson decays is the parametrization 
of the final state interaction effects by means of a resonant rescattering of the final 
particles, for B decays final state interaction effects are in general expected to be 
small. This seems rather plausible because the final decay products are moving away 
from each other with large momentum; due to the relativistic time dilatation, the 
formation time of the final particles is large and, when formed, they are far away from 
the color sources, which implies that strong phases induced by the color interactions 
should be small (see also ||). Some evidence for this comes for example from 
exclusive non-leptonic two-body B decays 0. 

The expectation that the rescattering effects in the final state due to soft inter- 
actions become negligible in the m# — ► oo limit has been challenged in ||. Taking 
B — ► mi and B — > Kix as examples, these authors assume Pomeron dominance and 
Regge theory [Tl| at high energy (i.e. large \fs = m B ) to estimate the size of the 
inelastic final state interaction and its effect on CP- violating asymmetries, such as 
Ar = T(B~ — > K~7t°) — T(B + — > K + 7r°); their conclusion is that these inelastic ef- 
fects, mainly induced by the Pomeron, are sizeable. Since a reliable way to compute 
these soft, non perturbative, effects is missing at the moment, the rather pessimistic 
conclusion reached in || is that the presence of final state interactions will limit 
the accuracy of the standard approaches to B decays that are based on the use of 
perturbative QCD supplemented by the factorization hypothesis, even though in the 
tub — > oo limit both methods are expected to be more and more reliable. 

A rather different conclusion has been reached, on the other hand, by Zheng in [|J . 
Using basically the same approach (i.e. an approximate evaluation of the final state 
strong interaction S"-matrix based on the Regge model) and considering B — > DK 
non leptonic decays, Zheng concludes that the approximation of neglecting final state 
strong phases is " accurate up to, roughly speaking, about 10%" . Because of these 
conflicting results, we believe worthwhile to investigate more accurately this problem. 
Therefore, in the present paper, we consider other decay channels: B — > Dp, D*n 
using the same model (Regge theory) already employed in J|] and ||. Apart from 
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the choice of a different decay mode, our treatment differs from || and because 
we explicitly compute some inelastic effects, i.e. the Dp — > D*tc rescattering. The 
two-body D*ir, Dp B decays are presently investigated by several experiments (for a 
review see and present a noticeable theoretical and experimental interest, related 
to the validity of the factorization approximation and to the sign of the ratio of Wilson 
coefficients 02/(21- 

For decay amplitudes, final state interactions are taken into account by means of 
the Watson's theorem f|, [L2| : 



A = ^SA b (f ) 

where S is the S'-matrix, A b are the bare amplitudes, i.e. decay amplitudes with 
no final state interactions, and A are the full amplitudes. The S"-matrix relates 
amplitudes with the same isospin I and a given total angular momentum J. Since 
we consider the amplitudes A[ = A(B — > D*tt)i and A 2 = A(B — ► Dp)?, the final 
state has J = 0; moreover the vector mesons (p or D*) in the final state can only 
have longitudinal helicity (A = 0), which means that, effectively, we are dealing with 
a situation analogous to a decay into scalar particles, which simplifies considerably 
the formalism. 

Let us begin by writing down explicitly the isospin amplitudes (/ = 3/2 and 
I = 1/2) in terms of the physical amplitudes 



A(B^D*n) 3/2 = -^/2/3A(B° ^D*\ ) + ^1/3A(B° ^D*-n + ) 

A(B -> D*ti) 1/2 = ^fl/3A(B° -> D*°7i°) + ^2/iA(B Q -> D*-ti + ) (2) 



and, similarly : 



A(B^Dp) 3/2 = -^2/3A(B°^D p ) + ^l/3A(B°^D-p + ) 
A{B^Dp) l/2 = Jl/3A(B° - D°p°) + J2/3A(B° - D~p+) . (3) 
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The J = ^-matrix, for each given isospin channel, must satisfy the unitarity 
relation. In standard notation, the two-body S'-matrix elements are given by = 
1,2) 

S% ]1 = Sa + 2iVw4 0)i > ( 4 ) 

where the J = 0, isospin / amplitude Afj = A« (s) is obtained by projecting the 
J = angular momentum out of the amplitude A^(s,t): 



Working in the approximation mp = mp*, m n ~ m p ~ 0, with yfs = tub, we 
have Xi = (s — m 2 D ) 2 , pi = \fK/ s ~ (s — m 2 D )/s, t- = and t + = — (s — m 2 D ) 2 / s. 

In order to compute Eq. (^) and Eq. (|5|), we need A^s, t). As we stressed already, 
in order to evaluate these amplitudes, we will work in the Regge model, which should 
be a reasonable theoretical framework due to the rather large value of s = m 2 B . In 
terms of the Pomeron (P), which is the leading contribution, and the non-leading 
trajectories p, /(1270) and 7r, neglecting Regge cuts Q, we have: 

Al{ 2 = A(D*Tr^D*7r) 3/2 = P + f + p 

A% 2 = Al{ 2 = A(D*n^Dp) 3/2 = n 

A% 2 = A(Dp^Dp) 3/2 = P' + f' + p' (6) 

and 
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A^ = A(D*n^D*n) 1/2 = P + f-2p 
A 1 / 2 = A 1 2 /2 = A(D*7r^Dp) 1/2 = -2<K 

Al /2 = A(Dp^Dp) 1/2 = P f + f-2p> . (7) 

We observe that the primed (P', /', p') contributions may differ from the un- 
primed ones only for a numerical coefficient; we also observe that the leading Regge 
trajectories in the off-diagonal matrix elements should be to and A 2 ; however they 
only contribute to the helicity-flip amplitudes (with A = ±1) that are not of interest 
here, which is why we take into account the next-to-leading trajectory, i.e. the pion 
(7r) Regge exchange. 

Let us first consider the Pomeron contribution, that we parametrize as follows 

P= - (3 p g(t){£) ap{ \-^ , (8) 

with so = 1 GeV 2 and 

a P (t) = 1.08 + 0.25t (t in GeV 2 ) , (9) 



as given by fits to hadron-hadron scattering total cross sections |R|, ||14] . The product 
f3 p - git) = (3 p {t) represents the Pomeron residue; for the t-dependence we assume 

[s — t/rrLpY 

which is motivated by the analogy with the electromagnetic form factor and by the 
smallness of t, due to the exponential damping in (s/ so) ap ^ . As for the residue at 
t — 0, i.e. j3 p , we assume, as usual, factorization: 

f = M (ii) 



1 A discussion on Regge cuts is contained in 



for the elastic D*tt — > D*ir amplitude. 

The residue at the vertex Ptttt can be extracted from proton-proton and pion- 
proton high energy scattering; we find 

* \f v = 5.1 (12) 

which is consistent with the hypothesis of the additive quark counting rule: (3 P = 
2(3 p (uu), (3 P = 3f3 p (uu). 

It is worthwhile to remark at this stage that we can obtain, from 7^ high energy 
scattering data, the Ppp residue (3 P by making the assumption of Vector Meson 
Dominance (VMD); in this way we find the approximate relation 

f P * ft (13) 



which is numerically valid within 15%. Eq. (|13"D is also consistent with the additive 
quark counting rule. As for the coupling of the Pomeron to charm, assuming again 
the additive quark model, one has 

ft = ft* = ft(cu) + (3 p (uu) (14) 

and for (3 p (cu) one has to assume as an input some theoretical ansatz; for example 
in H it is assumed: 

(3 p (cu) = ^P P (uu) . (15) 

We shall assume Eq. fll~5|) as well and will comment on this choice below. Let us 
observe that, by this assumption, we find 

P = P' , (16) 

a result to be used in Eq. (Q)-Eq. (0). 

Let us now consider the non-leading Regge trajectories R{= p,f,n). For these 
exchanges we write the general formula 

1 _|_ ( \s R -iira R (t) 

R = -(3 R { V T ( l R ~ «i ? (t))(a / ) 1 "^(«'^) aB(t) (17) 



as suggested in fl5| . an(t) is the Regge trajectory given by 

&R(t) = ocr(0) + a't , (1 

with an universal slope a' = 0.93 GeV~ 2 and a p (0) = a/(0) = 0.44 and a^(0) 
— a'm^. Ir is the lowest spin occurring in the exchange degenerate trajectory (l p 
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If = 1, l n = 0) and sr is the spin of the exchanged meson in the Regge amplitudes 
(s p = 1, S/ = 2, s n = 0). The choice Eq. fll7D for the Reggeized amplitudes is 
suggested by the high energy limit of a Veneziano amplitude. Since or(£) = sr + 
a'(t — m R ), near t = m R , Eq. (pT7|) reduces to 

^ ~ ^ r ^ (is) 

(mjj - t) 

which allows us to identify (3 R as the product of two on-shell coupling constants (as 
l"5f , we neglect here the t-dependence of the Regge residues^]). 



in 



To compute the different residues, we assume factorization, exchange degeneracy 
and SU (4) symmetry (which is of course largely violated, but should at least provide 
us with an order of magnitude estimate). Therefore we have [3 P = ^^b*d*i an d also 
Pdd = Ann an d Pp P — Pd*d*i Using as an input the g p7rn coupling constant and Vector 
Meson Dominance to relate j3 p pp to the electromagnetic coupling of the p + particle, we 
finally obtain j3 p = ~ 23 and p = p' , f = f ■ As for the pion exchange, we obtain, 
by this method, (3^ po = 4.6 GeV from p — > nn decay, and /3p D * = 1.7 GeV from a 
theoretical estimate of the D*Dir coupling ||16|| , with the result /3 n = 7.8 GeV 2 . Using 
these results, together with Eq. (f|)-Eq. (0), we obtain the following results for the 
J = 0, isospin 1 = 3/2 and 1/2 2x2 S'-matrices: 

^(0)3/2 = ( 0.76 - 0.06* -(0.17 + 1.6*) x lO" 2 \ 

V —(0.17 + 1.6*') x 10~ 2 0.76 - 0.06*' J 1 ' 

and 

^(0)1/2 = ( 0.71 - 0.02* +(0.34 + 3.2*) x 10~ 2 \ 

V+(0.34 + 3.2i) x 10~ 2 0.71 - 0.02*- ) 1 ' 

Let us now comment on these results. First of all, as one can see, S^ 3 ^ 2 and 
5'(o)i/2 j n fl2~0|) and Eq. ( pi]) are not unitary matrices, the reason being that 
other inelastic effects, besides the Dp — > D*n final state interaction, are present. We 
shall comment on this point later on; for the time being we observe that, because 
of the smallness of the off-diagonal matrix elements in Eq. (BO) and Eq. (2T), in- 
elastic Dp — > D*tt scattering is not expected to play a major role in determining, 
through Watson's theoremf] , the full (B — > f) amplitude. For example, the bare 
amplitude Ab(B° — > D°p°) does contribute, through final state interactions, to the 

2 We have verified that this does not alter our numerical conclusions. 
3 Wc find 

(VStfW 9 ) = (v 7 ^ 5 ^)^ = 0.87-0.03^; (v 7 ^ 3 / 2 )^ = (v 7 ^^ 3 ^)^ = -(0.07 + 0.9*) x lCT 2 
(VSW 1 / 2 } = (yW^j = 0.84-0.02* ; (VS^Va) = (^ s(0)1/2 ) 21 = +(0-16 + 1.9*) x 10" 2 
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decay amplitude A(B° — > D* ti ), but its contribution to the width does not exceed a 
few percent. Moreover, the bare amplitudes Ab(B° — > D*~n + ) and Ab(B° — > D~p + ), 
that, in principle, also contribute to A(B° — ► D*°7r°) via final state interactions and 
can destroy the simple predictions based on factorizationF], give a negligibly small 
contribution (of the order 2%) to the width. 

As mentioned earlier, the above results are obtained in the Regge model with ex- 
change degeneracy. One could avoid assuming exchange degeneracy and take different 
residues and intercepts for the u and / as trajectories fitted by total cross sections 
and elastic scattering data [13. In this way, one would obtain: = 9 p f = 21.5 
a p (0) = 0.57, a/(0) = 0.43, a' = 0.93 GeV~ 2 . The difference in the residues {(3 P ^ (3 } ) 
is largely compensated by the different intercept and the final results, as expressed 
by Eq. ( |20| ) and Eq. (^1|), would be unaffected (changes would be less than 3%). 

These results should hold, at least qualitatively, also if we enlarge the basis of 
eigenstates to enforce the unitarity of the S-matrix. To show it explicitly implies a 
considerable amount of work, but we can be convinced of it by going to the approxi- 
mation of neglecting the pion exchange contribution 7r in Eq. (|J) and Eq. fl7|) because 
of its smallness, and introducing, similarly to f|] and ||, two effective states |3 > and 
|4 > to take into account the inelastic scatterings of D*ir and D*p states respectively. 
In this way the unitary 4x4 S'-matrix can be constructed 



/ r]e 2tS iy/1 - ^ e ^ 5+5 ^ \ 

(0)7 _ r]e 2iS - ri 2 e l ( s+s ^ 

iy/1 - r]' 2 e i( - s+Sl) r]e 2lSl 

V iy/1 - n' 2 e i{5+52) r]e 2i52 I 

where, for / = 3/2, 7] = 0.76 and 6 = —0.04, while for I = 1/2, rj = 0.71 and 
5 = —0.02 (5 1 and 62 are free parameters in this model). For reasonable values of 5\, 
82 (Si ~ 5), we basically find the same result as in 0, i.e. that, within an uncertainty 
of 10 — 20%, the results that are obtained by including final state interactions do not 
differ from those obtained by the bare amplitudes with no final state interaction at 
all. 

We can apply these results to the calculation of the width of some two-body 
B decays into charmed states. With inelastic effects given mostly by the Pomeron 
contribution and neglecting the small strong phases S, the full decay amplitudes A 

4 In the factorization approximation they are proportional to the Wilson coefficient 01, while the 
bare amplitudes Ab(B° — > D°p°) and Ab(B° — > D*°tt°) depend on the much smaller coefficient 02, 
see e.g.0. 
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for B — > Z)*7r and -B — > Dp, according to Eq. ([!]), are (see also H): 



A = xA b + Vl - x 2 A[, (22) 



where is the B decay amplitude into an inelastic final state and x = y (1 + 77) /2. 
The decay rates with inelastic effects included can be obtained from 

\A\ 2 = \A b \ 2 + (l-x 2 )(\A' b \ 2 -\A b \ 2 ) . (23) 

Now, if A' b is comparable to A b , the decay rate thus obtained would be close to the 
rate obtained without final state interaction effects. This agrees well with experiment 



II |, [11] for neutral B decays into charged final states (e.g. D~ir + } D~p + , D*~n + ) 
where the factorization model works rather well for those amplitudes which do not 
receive contributions from the short-distance operator O2 which depends on a^- For 
those amplitudes which depend on both the short-distance operators 0\ and 2 , such 
as B + decays into D°ir + , D°p + , D*°tt + , the rates obtained in the factorization model 
without final state interaction effects are significantly smaller than the experimental 
results using the value a 2 — 0.11 as predicted by QCD. This suggests that in B decays 
there might be non-factorization contributions [|19|, |2(| . 

Which conclusions can we draw from our analysis? As a preliminary remark, let 
us observe that, even though we have dealt here with specific decay channels, some 
aspects of our analysis are general enough to be applied also to other decay modes. For 
example we have found that the inelastic amplitude D*n — ► Dp is strongly suppressed 
as compared to the elastic ones. This result does not depend on the dominance, in 
this case, of the pion exchange which has a small intercept, i.e. a T (0) ~ 0, but would 
hold quite generally because of the suppression of the non- leading Regge trajectories 
as compared to the Pomeron contribution. Therefore the only surviving inelastic 
effects at y/s = should consist of inelastic multiparticle production and should 
be dominated by Pomeron exchange. We have found, in agreement with ||, that 
also these inelastic effects should not destroy the predictions for exclusive two-body 
non-leptonic B decays based on factorization and perturbative QCD. In ||, charmless 
final states were considered and the opposite conclusion was reached. The apparent 
contrast has its origin in the hypothesis contained in Eq. (|T5|), that strongly reduces 
the Pomeron coupling to charm and therefore reduces the inelastic effects in all the 
channels with charmed mesons in the final state. This explains why [[5] finds significant 
final strong interaction effects in the considered channels. Eq. ( |15D has been assumed 
in § on a pure theoretical ground, since there is no experimental information on 
f3 p (cu). However a dependence on the quark mass can be clearly seen in the Kp 
total cross section which is asymptotically smaller than in np scattering ; a fit to the 



8 



asymptotic cross section is obtained by taking a reduction of 2/3 in the Pomeron- 
quark residue: (3 p (su) ~ |/3 p (wu); this shows a rather strong decrease with the quark 
mass. A powerlike dependence on the quark mass is certainly compatible with the 
result of Eq. fll5|), and, therefore, even if a clear numerical extrapolation is difficult, 
the assumption given by Eq. (|T5| ) is reasonable. Therefore we are led to the conclusion 
that final state interactions produce only moderate effects in the amplitudes provided 
that there are charmed (with or without open charm) particles among the decay 
products. On the other hand, for light charmless particles in the final state, the 
Pomeron contribution increases and rescattering effects may alter significantly the 
simple predictions based on perturbative QCD and heavy quark effective theory. 
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